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Abstract 

The isosinglet combination of the chiral-odd twist-3 distribution function e"(x) -|- e'^(x) of the 
nucleon has outstanding properties that its first moment is proportional to the well-known ttN 
sigma-term and that it contains a (5- function singularity at x = 0. These two features are insepa- 
rably connected in that the above sum rule would be violated, if there is no such a singularity in 
e"(x) -|- e'^{x). Very recently, we found that the physical origin of this 5- function singularity can be 
traced back to the long-range quark-quark correlation of scalar type, which signals the spontaneous 
chiral symmetry breaking of the QCD vacuum. The main purpose of the present paper is to give 
complete theoretical predictions for the chiral-odd twist-3 distribution function e"(x) of each flavor 
a on the basis of the chiral quark soliton model, without recourse to the derivative expansion type 
approximation. These theoretical predictions are then compared with the empirical information 
extracted from the CLAS data of the semi-inclusive DIS processes by assuming the Collins mech- 
anism only. A good agreement with the CLAS data is indicative of a sizable violation of the vrA'^ 
sigma-term sum rule, or equivalently, the existence of a 5- function singularity in e"(x) + e'^{x). 
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I. INTRODUCTION 



It is a widely accepted common belief now that the nonperturbative dynamics of 
QCD (chiral dynamics) is an indispensable element to understand high-energy deep inelastic 
scattering observables. Undoubtedly, the reconfirmation of this natural fact strongly owes to 
the two remarkable experimental discoveries in this field [l]-[3]. They are the unexpectedly 
small quark spin fraction of the nucleon revealed by the EMC measurement [1],[2] and 
the light- flavor sea-quark asymmetry confirmed by the NMC measurement [3]. The most 
successful theoretical studies of parton distribution functions have been carried out within 
the framework of the chiral quark sohton model (CQSM) [4] -[15], which is an effective 
model of baryons maximally incorporating the spontaneous chiral symmetry breaking of the 
QCD vacuum. In fact, we claim that it is so far the only one effective model of baryons 
which are able to explain the above two remarkable findings simultaneously within the single 
theoretical framework [16]-[19]. 

Very recently, we became aware of another novel example in which nonperturbative QCD 
dynamics play an unprecedented role in the physics of parton distribution functions. It 
concerns the possible existence of a delta-function singularity at the Bjorken variable x — 
in the chiral-odd twist-3 distribution function e{x) of the nucleon [20], [21]. This distribution 
function itself, together with its first moment sum rule giving the familiar ttN sigma-term, 
have been known for a long time [22]. In spite of several interesting theoretical features, 
however, this distribution function has been thought of as an academic object of study, since, 
because of its chiral-odd nature, it does not appear in the cross section formula of inclusive 
DIS scatterings. The situation changed drastically, however, since the CLAS Collaboration 
was able to get the first experimental information on this interesting quantity through the 
measurement of the azimuthal asymmetry Am in the elcctroproduction of pions from deeply 
inelastic scattering of longitudinally polarized electrons off unpolarized protons [23]-[25]. 

Some years ago, within the framework of perturbative QCD, Burkardt and Koike noticed 
that the first moment sum rule (or the ttN sigma-term sum rule) for e{x) holds only when 
e{x) has a 5-function type singularity at the Bjorken variable x — [26]. Unfortunately, 
the physical origin of this singular term is not very clear in this perturbative analysis. Very 
recently, two independent proofs were given to the fact that the physical origin of this 
5-function singularity can be traced back to the nonvanishing vacuum quark condensate 
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which signals the spontaneous chiral symmetry breaking of the QCD vacuum [20], [21]. An 
interesting question is whether we can verify experimentally the existence of this 5-function 
singularity in e{x). Unfortunately, the point x = is experimentally inaccessible. This 
means that, if there really exists such a 6{x)-tjpe singularity in e{x), the experimental 
measurement would rather confirm violation of this ttN sigma-term sum rule. Nonetheless, 
since e{x) in the region x ^ can in principle be measured, theorists are challenged to 
explain its behavior. 

The first theoretical study of e{x) was done by using the MIT bag model [27] . (See also 
[28].) However, this estimate based on the bag model cannot be taken as a realistic one 
by the following reasons. First, its prediction for the magnitude of the nN sigma-term is 
far from reliable. Second, more seriously, it cannot reproduce the (5-function singularity 
of e{x). Both these features (they are not actually unrelated) are easily anticipated, since 
the MIT bag model is essentially a relativistic quark model with Nc{— 3) valence quark 
degrees of freedom only, and the reproduction of the nonzero vacuum quark condensate is 
beyond the range of applicability of this model. The first realistic investigation of e{x) was 
carried out by Efremov et al. on the basis of the chiral quark soliton model but within 
the "valence" quark only approximations [29], [30]. More recently, the present authors and 
Schweitzer independently carried out more careful analysis of the contribution of the Dirac 
sea quarks on the basis of the gradient expansion type approximation and confirmed that the 
isosinglet combination of e{x) certainly contains 5-function type singularity [20], [21]. After 
some analysis of higher derivative terms of the expansion, however, Schweitzer retreated 
to the assumption that the contribution of the Dirac sea quarks is saturated by this 5{x) 
term alone. As admitted by himself, however, whether this last assumption is justified 
or not is far from trivial [20]. To confirm it, one has to carry out an exact numerical 
calculation within the model without recourse to the gradient expansion type approximation. 
Furthermore, to compare the predictions of the model with the experimental data of the 
CLAS collaboration, one must know e"(x) of each fiavor a. To this end, only the knowledge 
of the isoscalar combination e"(a;) + e'^{x) is not enough. We need another independent 
combination, i.e. the isovector distribution e'^{x) — e'^(x). Within the framework of the 
CQSM, this latter distribution survives at the next-to-leading order in l/N^ expansion and 
it was left untouched in [20]. 

In view of the above-mentioned circumstances, we think it important to carry out an exact 
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model calculation within the CQSM for both of the isoscalar and isovector combinations of 
the chiral-odd twist-3 distribution function e{x). We also think it useful to analyze the first 
and the second moment sum rule for e"(x) +e''(x) and e^{x)—e'^{x) within the CQSM in the 
light of the corresponding sum rule expected in the general framework of perturbative QCD. 
The predictions of the model for e"(a;) and e'^{x) (as well as the corresponding distributions 
for antiquarks) are then used as initial distributions given at the model energy scale around 
600 MeV (or Q"^ ~ 0.30 GeV^), and they are evolved to higher Q'^ for the sake of comparison 
with the phenomenological information obtained by using the CLAS measurement. 

The paper is organized as follows. In sect. 2, after a brief introduction of the basic idea 
of the CQSM, the theoretical expression for e"(a;) + e'^(x) and e"(j;) — e''{x) arc given. The 
fundamental moment sum rules for these distributions arc also discussed here in some detail. 
Sect. 3 is devoted to the discussion of the numerical results. Finally, in sect. 4, we summarize 
what we have found. 



II. e{x) IN THE CHIRAL QUARK SOLITON MODEL 

The chiral-odd twist-3 quark distribution e"(a;) of flavor a inside a nucleon with 4- 
momentum P, averaged over its spin, is defined by 

/oo fj y~ , 

where -0^ ^re quark fields. Similarly, the corresponding antiquark distribution is defined as 

/CXD f\ 'y~ 
-oo ZTT 

with -0*^ being the charge-conjugate field of '0- Here we use the standard light-cone coordi- 
nates 

The variable x denotes the Bjorken variable, x = —(fj {2P-q), with q being the 4-momentum 
transfer to the nucleon. Taking account of the charge-conjugation property of the relevant 
quark bilinear operator, one can formally extend the domain of quark distribution functions 
from the interval 0<j;<l to — l<x<l, such that 

e^(x) =e"(-x), (0<x<l), (4) 



4 



which dictates that the distribution function with negative x should be interpreted as anti- 
quark one. 

Although the above definitions of the quark and antiquark distribution functions are 

frame-independent, it is convenient to perform the actual calculation in the nuclcon rest 
frame. In this frame, we have P"'" = Mjv/\/2, and the distribution function is reduced to 

Throughout the paper, we will confine ourselves to two flavor case of u- and d-quarks, 
and neglect strangeness degrees of freedom in the nucleon. Consequently, we have two 
independent distributions, i.e. the isosinglet distribution e^'^^^\x) = e"(a;) -|- e^{x) and 
the isovector one e^'^^^^x) = e"(x) — e'^{x). In the case of e^'^^^^x), we simply sum up 
(5) over the flavor components. On the other hand, for e^^^^\x), we have to sum up the 
representation after inserting T3 matrix in (5). 

For obtaining quark distribution functions, we must generally evaluate nucleon matrix 
elements of bifocal and bilinear quark operators containing two space-time coordinates with 
light-cone separation. The startingpoint of our theoretical analysis is the following path in- 
tegral representation of the matrix elements of a bifocal and bilinear quark operator between 
the nucleon states with definite momentum P : 

{N{P)\iP\0)^''iP{z)\N{P)) ^ ^ J d^x d^y e-'^ '' e'^ y J VU 



xe-x.Y>[i j d'^x ip {i^ - MW)ip], (6) 



where 



C = ^ - MW'{x))'il), (7) 

with 

C/^n^)=exp[^75r-7r(a;)/A], (8) 
being the basic lagrangian of the CQSM with two flavors. The quantity 



5 



is a composite operator carrying quantum numbers J J3, TT3 (spin, isospin) of the baryon, 
where ai the color indices, while ^jX^tts ^ symmetric matrix in spin flavor indices /j. 
We start with a stationary pion field configuration of hedgehog shape : 

7r{x)^UrF{r). (10) 

Next we carry out the path integral over 7r{x) in a saddle point approximation by taking care 
of two zero-energy modes, i.e. the "translational zero-modes" and "rotational zero-modes" . 
Under the assumption of "slow rotation" as compared with the intrinsic quark motion, 
the answers can be obtained in a perturbative series in Q, which can also be regarded as 
a l/Nc expansion. Up to the first order in the collective rotational velocity Q, the only 
surviving contribution to e^'^^'^\x) arises at the term of this expansion, since the 

0{Q^) term vanishes identically due to the hedgehog symmetry. On the other hand, the 
first nonvanishing contribution to e'^'^^^^x) arises at the 0{^l^), since the leading 0{^l^) 
contribution vanishes due to the hedgehog symmetry. Then, between the magnitude of the 
above two distributions, one may expect the following large- A^c relation : 

|e"(a;) + e'^{x)\ - |e"(x) - e'^{x)\ . (11) 
A. Isosinglet distribution e(^~'')(£c) 

The isosinglet combination of the chiral-odd twist-3 unpolarized distribution is given by 

eM(a;) ^ e^{x) + e''{x) = T ^ e'^''-^^ {N\mm\N)L=-.o,z.=o. (12) 

J— 00 ZTT 

Following the general formalism developed in [4], [5], [9], the isosinglet distribution in the 
CQSM is given in the following form : 

e(^=o)(x) = -TVcMjv Yl {nh^SixM^-ps - En)\n) (13) 

n>0 

= iV,M^^(n|705(a;M^-p3-^n)|n), (14) 

n<0 

where, \n) and En are the eigenstates and the associated eigenenergies of the static Dirac 
Hamiltonian 

H^-icx-V + (5Me'^^'^-^''^'\ (15) 

with the hedgehog background. Here, the summation X^„<o in (14) is meant to be taken 
over the valence-quark orbital (it is the lowest energy eigenstate that emerges from the 
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positive-energy Dirac continuum) plus all the negative-energy Dirac-sea orbitals. On the 
other hand, the summation J2n>o ™ (13) is meant to be taken over all the positive-energy 
Dirac continuum excluding the discrete valence orbital. We recall that the CQSM is defined 
with some appropriate rcgularization. In fact, without regularization, e^^^^\x) is quadrat- 
ically divergent, and no practical meaning can be given to either of (13) and (14). The 
ideal regularization scheme for our purpose is the Pauli-Villars subtraction scheme, since it 
preserves several fundamental conservation laws of field theory [4], [5]. Furthermore, it is also 
expected to preserve the equivalence of the two ways of computing the quantity in question, 
by using (13) and (14). In the present study, we use the double subtraction Pauh-Villars 
scheme as introduced in [31], since e^'^^'^^x) divergees like the vacuum quark condensate. 
In this scheme the distribution e*^^""^ (x) is replaced with a regularized one defined as 



e 



(^=°)(a;) ^ e(^=°)(a;)^ - [^) e^^='\x)^'. (16) 

Here e(x)^* is obtained from e(x)'^ by replacing the mass parameter M by Aj. It was shown 
in [31] that, if the parameters ci, C2, Ai, and A2 are chosen to satisfy the two conditions : 

2 /A\2 



1=1 



l-gc.(^)*^0, (18) 

the quadratic as well as the logarithmic divergences in the vacuum quark condensate are 
completely eliminated. 

Actually, we are interested in the nucleon observables measured in reference to the phys- 
ical vacuum, so that e*^^"°^ (x) should be replaced by 

e(^=o)(^) ^ e(^=°)(x) = e^J='\x) - e^u=i\^)- (19) 

Here the vacuum subtraction term e^'~^\x) is obtained from e^'~^\x) by setting U — 1 
or F{r) = 0, and by excluding the sum over the discrete valence level. We point out that, 
due to the energy-momentum conservation embedded in the factor ^(xMjv — Ps — En), the 
vacuum subtraction terms are required only for x < in the occupied form (14), and for 
X > in the non-occupied form (13). This means that the vacuum subtraction terms need 
not be considered when e^-^^^^x) is evaluated in the following way, i.e. if it is evaluated by 
using the occupied form for x>0, while by using the nonoccupied form for x<0. 



1. momentum sum rules of e^^ 



The most important information of the distribution functions are generally contained in 
their first few moments of lowest orders. This is also the case for the distribution e^^"'^) (x) . 
In a recent paper, Efremov and Schweitzer reviewed some of the important sum rules for the 
chiral-odd twist-3 distribution functions in an enlightening way [32]. Their argument starts 
with the general definition of the distribution with flavor a as 

''"^'^^ = 2i^/ ^^'''Wa{m,>^ri]i:a{Xn)\N), (20) 

where [0, An] denotes the gauge hnk. By using an operator identity following from the QCD 
equation of motion, e^{x) is shown to be decomposed in a gauge invariant way into the 
following three pieces : 

e"(^) = e%Jx) + et^,{x) + e^_(x). (21) 
Here, e^j„^(a;) denotes a singular term given by 

e:,^^{x) = 5{x){N\rr\N). (22) 

On the other hand, e^^.^{x) is a genuine twist-3 part of e"(x) that contains information on 
quark-gluon-quark correlations. Finally, e'^nassi^) denotes the term arising from the nonzero 
current quark mass. It is a somewhat peculiar function defined through its Mellin moments 

as [33]- [36] 

£ x'^-^ el,,,Ax) dx = 5„>i • ^ £ x^-' mx) dx, (23) 

with fi{x) being the twist-2 unpolarized distribution with flavor a. The presence of the 
factor Sn>i here dictates that the flrst moment of e^„^^(x) vanishes, 

£^C,ssi^)dx^O. (24) 

It is also known [33]-[36] that the first two basic Mellin moments of ^^^^{x) vanish, i.e. 

£^ x"~^e^^3(a;)dx = for n = 1, 2. (25) 

Putting the above-mentioned properties altogether, the first moment sum rule for the 
isoscalar combination of e"(x), i.e. e*^^"°^(x) takes the form. 

eM(x)da; = ^, (26) 
J-i mo 



which is nothing but the nN sigma-term sum rule. Note that this sum rule is saturated by 
the first term of (21) alone. On the other hand, the second Mellin moment of e^'^^°^(x) is 
given by 

j\e^^='\x)dx = ^ ■ N,, (27) 

where Nc is the number of color, which coincides with the number of quarks contained in a 
baryon-number-one system, i.e. Nc — 3. We point out that this second Mellin moment of 
e(^=°) vanishes in the chiral limit of mo — 0. 

Next, we turn to the discussion of the moment sum rule in the CQSM. Integrating (14) 
over the first moment of e*^"^^^'' (x) is given as 

f e^^=^\x)dx = iVe ^(n|7». (28) 

n<0 

Since the r.h.s. of this equation is nothing but the scalar charge a of the nucleon within the 
CQSM, the sigma-term sum rule immediately follows 



e^'^=^\x)dx^a^^. (29) 
mo 



The way of this sum rule being satisfied is far more delicate in the CQSM than in the above 
QCD-motivated analysis. As shown by our previous study, although the model certainly 
predicts the 5(a;)-type singularity in e^'^^^\x), this term alone docs not saturate the ttN 
sigma-term sum rule. The model also predicts nontrivial structure of e^'^^'^^x) at a; 7^ 0, 
which may contribute to the first moment sum rule. We shall discuss this point in more 
detail in the next section. 

Turning to the second moment, it is easy to show from (14) that 



/_ 



xe^'^='\x)dx = ^ E("l7°(P3 + E^)\n). (30) 



N „<0 



Owing to the hedgehog symmetry of the soliton, the term containing ^y^ps vanishes, and we 
are left with 



xe 



(r=o). 



N, 



^)dx=-^T.En{n\7». (31) 



Following [20], it is convenient to rewrite the r.h.s. of the above equation in the following 
manner. First, notice the identity 

K(n|7» = l{n\{H,^'}\n) = mo + M{n\^{U + C/t)|n). (32) 



Here, we have tentatively restored the current quark mass term in the model Hamiltonian 
H, just for the sake of explanation here only, i.e. we have used here 

H=-icx-V + /?Me^^5T.fF(r) ^ (33) 



Then, the second moment sum rule in the CQSM takes the following form 
with 



[' xe^^=''\x)dx = ^ (mo + (3M), (34) 

7-1 Mat 



f3^Y:{n\l{U + U^)\n). (35) 

n<0 ^ 

It is clear now that r.h.s. of this sum rule does not vanish even in the chiral limit of mo = 0, 
in contrary to the sum rule derived from the QCD-equation-of-motion method. We shall 
return to this point in the next section. 



B. Isovector distribution e^^~^\x) 

The isovector distribution is defined by 

e^^='\x) ^ e-{x) - e'{x) = H ^ e"-^^-^»(iV|^(0)r3^(z)|iV)|,3=_,„,,^=o. (36) 

Within the framework of the CQSM, e^'^^^\x) survives only in the next-to-leading order 
in the collective angular velocity Q. Following the formalism derived in [9], [10], the final 
answer is written in the form : 

a=l m=all,n>0 \ m n N / 

= {2n),M^^\ j: (n|r„|m)(m|r„/ f - ^ S'\ \n), (37) 

a=l m=all,n<0 \hm - i^n ) 

with (5„ = b{xM^ - En - p^) and S'^ = -§-^5{xMn -En- p^). Here / in the r.h.s. of (37) is 
the moment of inertia of the soliton, given by 

J -/Vc V- V- {n\Ta\m){m\Ta\n) 

I^-^l^l^l^ -^-^ ■ ^38) 

" a=lm>On<0 "^"^ 

In (37), {0)p should be understood as an abbreviated notation of the matrix element of a 
collective operator O between the (spin-up) proton state, i.e. 

= /*T=r3=l/2;J=J3=l/2MO[eA]*T=T3=l/2;J=J3=l/2MrfeA 

= (p,53 = 1/2101^,-53 = 1/2). (39) 
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In the present case, we have (ST'a)^ = 1. 

We immediately notice that the above expressions are not suitable for the actual numerical 
calculation. Here, we shall proceed as in the previous studies [9], [10]. First, note that the 
term containing the x-derivative of the 5-function in (37) can be rewritten as 

^.(-) = E E (n|..h)(.|../^„|„> ^ (40) 

o m=all,n<0 " 

Here, we have made of the completeness of the eigenstates \n) of the static Dirac Hamilto- 
nian H. (We recall that 62 (x) term originates from the non-locality in time of the operator 
'ip{Qi)T°"<p{z) in (36).) It should be recognized that the x-derivative of the isosinglet distri- 
bution e^'^^'^\x) appears in the right-hand side. Since we already know that the isosinglet 
distribution e*^"^^"-* (a;) has the 5{x) type singularity connected with the nonvanishing vacuum 
expectation, it then follows that 62 (x) has the derivative-of-5(x) type singularity. However, 
it is unlikely that the net isovector distribution e^'^^^\x) has such a singularity, because 
the QCD vacuum should not violate isospin symmetry so that vacuum quark condensate of 
isovector type must simply vanish. This apparent discrepancy can be resolved as follows. 
We first divide the double sum of (37) into the sum over terms with Em = and with 
Ejji 7^ En. The point is that the sum with E^, = E^ in ei{x) can be rewritten in a similar 
form as the corresponding term in 62(2;), 

ei(x) = Mn^^Y. J2 ^ {n\Ta\m){m\Taj°6n\n) 

a m=all,n<0 -^m ^71 

J2 {n\Ta\m){m\Ta'y'^dn\n), (41) 
ax o Q, ^<o,n<o 

(Em=En) 

e2{x) = -^-TT^IZ {n\ra\m){m\Ta-f^5n\n) 

ax 4i O „ m=all,n<0 

E {nK\m){m\Ta^'W- (42) 

ax ^1 O (J m<0,n<0 

Now, just as argued in [10], [9], E^. — E^ contribution in the double sums in ei{x) and 62 (x) 
precisely cancel each other. After regrouping the terms in such a way that this cancellation 

occurs at the level of analytical expressions, the 0{Vt^) contribution to the distribution 
function e*^^^^'*(a;) = e^{x) — e'''{x) can finally be written in the following form: 



e^'^=^Hx) 
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The fact is that, in the double sum of (42), the singularity connected with the non-zero 
vacuum quark condensate comes only from E^, — En contribution, i.e. the second term of 
(42). As mentioned above, after the E^ — E^ contributions in ei{x) and 62(0;) are canceled, 
these singularities disappear in (43). The final theoretical formula (43) is therefore free from 
any singularity which contradicts the symmetries of the QCD vacuum, and it provides us 
with a sound basis for numerical calculation. 

1. first moment sum rule of e^^^^\x) 

Here we discuss the first moment sum rule of the isovector distribution. Integrating (36) 
over X, we obtain 

J^^ e^^=^\x)dx = f^ie'^ix) - e^{x))dx = {N\i^T^i^\N) . (44) 

(Here, '0T3'0 should be taken as an abbreviated notation of / ip{y)T3ip{y) d^y, which gives 
the isovector scalar charge operator.) An interesting observation is that the first moment 
of e^^^^^ (x) is related to the non-electromagnetic mass difference of neutron and proton. In 
fact, the nonelectromagnetic neutron- proton mass difference is thought to be generated by 
the isospin breaking term in the QCD Hamiltonian : 

AH = Hhi^ (^^^^ _ ^,^,) . (45) 

Because of the smallness of all the masses m„, ma, md—mu compared with the typical energy- 
scale of hadron physics, we can treat Aif as the first order perturbation, thereby being led 
to the following formula for the non-electromagnetic mass difference between neutron and 
proton : 

{M„-Mp)QCD = {n\AH\n) - {p\AH\p) 

= {nid- mu){p\i>uiJu- ^/Jd^l^dlp): (46) 

where use has been made of the isospin symmetry for the unperturbative state \p), \n) (i.e., 
the invariance under the interchanges p n and u d). Empirically, the neutron- proton 
mass difference of QCD origin can be estimated from the observed mass difference by taking 
account of the correction due to the electromagnetic interactions : 

(M„ - Mp)QCD = {Mn - Mp)e^p - {Mn - Mp)e.m. ■ (47) 
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Using the values (M„ - Mp)exp ^ 1.29 MeV, (M„ - Mp)e.m. ^ (-0.76 ± 0.30)MeV [37], we 
obtain 

(M, - Mp)QCD ^ (2.05 ± 0.30) MeV , (48) 

To extract the first moment of e^'^^^^x) empirically, we need to know the value of — m^. 
By using — ~ 5 MeV, as an order of magnitude estimate, we obtain 

e(^=^) (x) dx = ~ ^^^^^^ 0.41 ± 0.06 . (49) 
J-i rricL - ruu 

On the other hand, the theoretical expression for the first moment of e*^^^^-* (x) is obtained 
from (43) as 



^ -J a n<Om>0 ^"i J^n 



Here, we have used the fact that, since the contribution e2{x) is a total derivative, it does 
not contribute to the integral of (50). After integration over x, the double sum over levels in 
(50) is naturally restricted to include only transitions from occupied to non-occupied states. 
This is reasonable, since the operator appearing in the r.h.s. of (50) is a local operator, and 
transitions from occupied to occupied states would violate the Pauli principle. Within the 
framework of the CQSM, we can evaluate the r.h.s. of (44), i.e. the isovector scalar charge 
of the nucleon {Nlipr^iplN), directly without passing through the distribution function. 
Since the resultant expression of {N\ipT^ip\N) precisely coincides with the r.h.s. of (50), we 
conclude that the first moment sum rule of e^-^"^^ (x) is properly satisfied within the model. 



III. NUMERICAL RESULTS AND DISCUSSION 

The numerical method used for evaluating e{x) in this paper is essentially the same as 
the one used for the computing the twist-2 distributions q{x), Aq{x), 6q{x) [8, 9]. The 
eigen-energies and the eigen-vectors of the static Dirac Hamiltonian H with the hedgehog 
background are obtained by diagonalizing it with the so-called Kahana-Ripka plane-wave 
basis [38] . Following them, the plane- wave states, introduced as a set of eigenstates of the 
free Hamiltonian Hq — —ioc • V-|- /3M, is discretized by imposing an appropriate boundary 
condition for the radial wave functions at the radius D chosen to be sufficiently larger 
than the soliton size. The basis is made finite by retaining only those states with the 
momentum k satisfying the condition k < kmax- As a results of using this discretized 
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momentum basis, a resultant distribution becomes a discontinuous function of x, due to the 
factor 5{xMn — En —ps)- In order to get a continuous function with a discretized basis, we 
introduce a smeared distribution function in the variable x as [5] 

1 (x-x')'^ 

e-y(x) = — -= / e e{x')dx', (51) 

■Jy/TT J-oo 

with a small but finite value of 7 (7 <^ 1). The smeared distribution is expected to be 
continuous when the separation between the discretized momenta is much smaller than the 
smearing width 7. Since the physical distribution corresponds to the limit 7^0, this forces 
us to employ a very large box size D to get a continuous distribution function. 

This procedure works very well at least for the standard distributions investigated so 
far. However, in the numerical calculation of e^'^^^\x), we have a new problem which we 
have not encountered before. Our expectation is that, if a S{x)-type singularity really exists 
in e^-^^^^x), the corresponding smeared distribution would have a Gaussian peak centered 
around x — with the width 7. The problem here is that the distribution function in 
question may also have a piece that is non-singular for all value of x. One might think 
that the contribution of the singular part can be disentangled from the total contribution 
by using the "unsmearing method" described in [5]. This is not feasible, however, by the 
following reasons. First, although the smearing procedure defined with (51) preserves the 
integral value of the distribution, we have no ad hoc way to know the overall coefficient of 
the 5{x) term of the distribution. Secondly, the small x behavior of the nonsingular part of 
the distribution would be hard to know, because it is buried in the very large contribution of 
smeared 6 function singularity. This point will be discussed in more detail in the following 
subsection. 

A. isosinglet distribution e(^~^)(£c) 

In the numerical calculation, we fix the pion weak decay constant /tt in (10) to its 
physical value, i.e., = 93MeV, so that only one parameter of the model is the dynamical 
quark mass M, which plays the role of the coupling constant between the pion and the 
effective quark fields. Through the present analysis, we use the value of M = 375 MeV, which 
is favored from the phenomenology of nucleon low energy observables. With M — 375 MeV, 
we have Ai ~ 627 MeV and A2 ~ 1589 MeV from the conditions (17), (18). The static sohton 
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energy obtained with these parameters is about 1018 MeV. We point out that, although the 
sohton mass emerges about 8% larger than the observed nucleon mass Mat, the consistency 
with the energy-momentum sum rule of the unpolarized distribution functions enforces us 
to use this value for Mjy in the following evaluation of the distribution functions. 

We start with showing the numerical equivalence of the final answers based on the non- 
occupied representation and the occupied one. The problem here is the dependence on the 
cut-off momentum kmax, which is introduced to make finite the discretized Kahana-Ripka 
basis set. Since the distribution e^'^^^\x) is ultra-violet finite after introduction of the 
double-subtraction Pauli-Villars regularization, one might expect that the answers would 
be stable as far as one takes k^ax much larger than the second Pauli-Villars cut-off mass 
A2 — 1.6 GeV. This is not the case, however. As clarified in [21], the 5-function type 
singularity in e*^^"°^(a;) is generated by the contribution of the infinitely deep Dirac-sea 
levels, which are naturally contained in either of the three terms, i.e. the main term and 
the two Pauli-Villars subtraction terms. This implies that the singularity, which will appear 
in the smeared distribution as a Gaussian peak around x — with the width 7, would be 
reproduced only in the ideal limit of kmax -^00. To achieve this ideal limit, we therefore 
use an extrapolation method explained below. For this extrapolation to be done smoothly, 
we first introduce an energy cut-off into the level sums (13) and (14) of the form. 

n>0 

[e"(x) + e<^(x)]^,^^,,, = N,M^ ^ (n|7°5(xM^ -f- E^)\n)R{Er,). (53) 

n<0 

Here, R{En) is a smooth regulator function with an energy cut-off, E^ax = \l k^ax + 
For this regulator function, we employ here a Gauusian function 

(54) 

following Diakonov et al. [5]. We first compute the level sums (52) and (53) for several 
values of kmax-, in the case of masses M, Ai, and A2 respectively, and then perform the Pauh- 
Villars subtraction, and finally remove the energy cut-off by the numerical extrapolation to 
infinity pointwise in x. In the present study, we use five data (corresponding to kmax/M = 
12, 16, 20, 24, and 28) and perform a least square fit of these data by using a fourth order 
function of 1/kmax- 



R{En) = exp -{En/E, 
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extrapolation 
l<™x/M = 28 
WM = 24 

KJM = 20 

i I k™x/M = 16 

■ ' k™x/M = 12 




FIG. 1: The kmax dependence of the Dirac-sea contribution to e^-^ based on the occupied 

representation. The sohd curve represents the extrapolated result. 

Now we are ready to show in Fig.l the kmax dependence of the Dirac-sea contributions 
based on the occupied representation for all values of x. Here we use a value of 7 = 0.1. This 
figure shows that the peak positions of the Gaussian-like function obtained with the finite 
cutoff energy deviate to the negative x region from the origin x = 0. This deviation of the 
peak position in the smeared distribution may be understood as follows. First, when one uses 
the occupied representation, the vacuum substraction as represented by (19) is necessary only 
for the region x < 0, while it is not necessary for 2; > 0, since the vacuum term identically 
vanishes for x > due to the restriction of the factor S{xMn — En —ps). Secondly, we recall 
the fact that the singular term of e^^"'^^ (x) emerges as a dehcate cancellation of two large 
numbers or the infinities, i.e. the difference between the main contribution with hedgehog 
background and the vacuum subtraction term obtained with U = 1. These two facts indicate 
that the use of the occupied form with some finite value of 7 can reproduce the redistribution 
of the delta-function strength at x = in the a; < region only, but it cannot do it properly 
in the x > region, as far as the finite energy cutoff is used. This is the reason why the 
Gaussian-like peak of the smeared distribution is shifted to the negative x region. One can 
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however confirm the behavior that the position of the Gaussian peak approaches a; = as 
the energy cutoff is increased. And, finally, with the extrapolation method, we obtain a 
reasonable result which shows that the peak of the smeared distribution is positioned just 
at a; = 0. (In the above analysis, we fix the box size to be DM = 20. As a matter of 
course, to get a physically acceptable answers, we must also investigate the dependence of 
the answers on the box size D. We found that, above DM — 20, the change of the small 
X behavior of e^'^^^\x) as illustrated in Fig.l is almost due to the increase of kmaxi ^-nd 
the answer is stable against the further increase of DM above 20.) After carrying out the 




X 

FIG. 2: Comparison of the Dirac-sea contributions to e^'^^^\x) based on the occupied (solid curve) 
and non-occupied (dashed curve) representations. 

similar analysis, this time, with use of the non-occupied representation, we can now compare 
the final numerical results for the Dirac-sea contribution obtained with the two alternative 
representations. Fig.2 shows this comparison. A reasonable agreement between the two ways 
of evaluating e'^'^^'^\x) confirms the equivalence of the two representations. At the same time, 
the analysis above establishes the existence of the 5-function singularity in e"(x) -|- e'^(x) on 
the numerical ground. Some difference between the two curves at the positive and negative 
X tails of the Gaussian like distributions would be a spurious one generated by the numerical 
extrapolation method. The contributions based on the occupied representation for x < 
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and the non-occupied representation for x > can be obtained after the cancellation of two 
large numbers, i.e., the main contribution and the corresponding vacuum subtraction term. 
On the other hand, if one uses the occupied representation for x > and the non-occupied 
representation for a; < 0, one is free from the spurious contribution due to the cancellation, 
so that the extrapolated curves at these tail region have reasonable smooth behavior. 

Although we were able to confirm the existence of S{x)-type singularity in the numerical 
analysis of e*^^=°^(a;), we cannot exclude the possibility that the e^'^^°\x) may also contain 
a regular term which is smooth in all the range of x. Is it possible to disentangle such 
non-singular term of e^'^^^\x) from the total contribution containing the singular one? One 
should recognize that it is not so easy by the following reasons. First, the dcconvolution 
method as proposed by Diakonov et al. does not work because of the very delicate nature 
of the singularity [5]. Second, we have no ad hoc way to know the coefficient of d{x) term in 
the original unsmeared distribution. Nevertheless, we found that the following trick works 
for obtaining the non-singular distribution excluding the 5{x) term. That is, as repeatedly 
emphasized, by using the non-occupied expression for x < and the occupied one for x > 0, 
we can avoid the vacuum subtraction. Interestingly, this also works to remove the singular 
contribution in the bare distribution, and the corresponding smeared distribution would not 
contain the Gaussian peak corresponding to S{x)-type singularity. (One should remember 
the fact that the vacuum term plays an indispensable role in reproducing the 5-function 
singularity in e^^"°^(x).) 

Fig. 3 shows the kmax dependence of the Dirac-sea contributions based on the occupied 
representation for x > and the non-occupied representation for a; < 0. One finds that 
the large and positive Gaussian peak, the reminiscence of the (5-function singularity in the 
bare distribution, does not appear any more. One can also see that the negative large 
contributions of the Dirac sea in the small x region tend to decrease as the cutoff momentum 
kmax increase. We again remove the energy cutoff by the numerical extrapolation to infinity 
pointwise in x. We observe some difference from the previous case, however. Owing to the 
feature that the 5-f unction singularity is already excluded in the present way of calculation, 
the kmax dependence is well reproduced by the linear function of l/k^ax as illustrated in 
Fig. 4. After this extrapolation procedure, the result shows a smooth behavior in the whole 
region of x except the region < 0.06 in which the answer is thought to contain some 
numerical instability generated by the extrapolation method. Neglecting the data in the 
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FIG. 3: The kmax dependence of the Dirac-sea contributions to e^^^°^(x) based on the occupied 
representation for x > and the non-occupied representation for x < 0. The sohd curves represents 
the extrapolated result. 

\x\ < 0.06 region, we make this extrapolated result smooth. After deconvoluting the smeared 
distribution with use of the Fourier and its inverse transforms, we obtain the final prediction 
for the distribution 6^-^"°^ {x) within the framework of the CQSM, the normalization point 
of which may be interpreted as about 600 MeV. 

Summarizing our analysis up to this point, the isosinglet part of the chiral-odd twist-3 
distribution is given as a sum of the valence quark and Dirac-sea quark contributions, 

e^^-'\x)^eZr\x) + eg:'\x), (55) 

where the Dirac-sea contribution consists of the singular term and the nonsingular (regular) 
term as 

el-°'W=CT(x)+eM(x), (56) 




I I extrapolated 

ll| KJM = 30 

f KJM = 28 

I KJM = 26 

^ WM = 24 



Shown in Fig. 5 are the final theoretical predictions for e^^^^^ (x) obtained in the above- 
explained way. The dashed curve here represents the contribution of Nc valence level quarks, 
while the dotted curve does the regular part of Dirac-sea contribution. The sum of these two 
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1 / k^x 

FIG. 4: The kmax dependence of es^a^\x) at a; = 0.12 and its linear extrapolation to kmax oo- 

contributions is shown by the sohd curve. (We recall that the 6{x)-tjpe singular term is not 
shown in this figure.) One can convince that the regular part of the Dirac-sea contribution 
shows a nontrivial structure in the x ^ region. 

After performing the numerical integration of the above distributions over x, one can 
obtain the contributions of the valence quark term and the regular part of the Dirac-sea 
term to the first moment sum rule : 

iyZr'\x)dx ^ 1.7, (57) 
j^^egf\x)dx ~ 0.18. (58) 

Note that the regular part of e^^a^\x) gives small but non-zero contribution to the sum 
rule. To determine the coefficient of the singular term in (56), we use the first moment sum 
rule (28) or (29) for e^^=°^(x), which was already shown to hold within the framework of 
the CQSM. We first recall that the r.h.s. of the sum rule (28) or (29) is the nucleon scalar 
charge defined by 

a=(iV|^„^„ + ^,^,|iV). (59) 
The point is that this low energy observable can be calculated within the CQSM, without 
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FIG. 5: The final theoretical predictions of the CQSM for e^'^^'^\x). The dot-dashed curve rep- 
resents the contribution of Nc valence level quarks, the dashed curve the non-singular part of the 
Dirac-sea contributions, and the solid curve does their sum. The 5-function singularity at x = in 
the Dirac-sea quarks part is not shown in this figure. 

asking for the distribution function e^'^^'^^x). It is given as 

^ = ^val + ^sea, (60) 



with 



Numerically, we find that 



so that 



^val = A^c(0|7», (61) 
asea = iVei:(n|7». (62) 

n<0 



CTval - 1-7, (Tsea - 10.1, (63) 



Cr = CTval + CTsea - H-S. (64) 

Then, by admitting the validity of the first moment sum rule, one can extract the coefficient 
of the 5 (x) term as follows : 

C = asea - e(^=°) {x)dx ~ 9.92. (65) 
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Our procedure for obtaining the coefficient C is different from tliat of Schweitzer [20] . After 
some consideration based on the gradient expansion analysis, he assumed that the Dirac-sea 
contribution to e^'^^^\x) is saturated by the 5{x) term with the coefficient E^rTv/'^^O; and 
simply neglected the possible existence of the nonsingular contribution. In his treatment, 
then, the nontrivial shape of e^'^^^\x) at ,x 7^ solely comes from the contribution of 
valence level quarks. Thus, the total distribution consists of these two terms as 

mo 

(Here, for simplicity, we ignore the term proportional to the product of mo and the unpolar- 
ized distribution function.) In our opinion, this procedure has a danger of double counting. 
Within the framework of the CQSM, the total vriV sigma term divided by the current quark 
mass mo is nothing but the total scalar charge a of the nucleon, which is made up of the 
two terms, a^ai and Osea- The x-integral of (66) would then lead to 

J ^ e^^=°\x)dx = {a^al + ^sea) + ^val, (67) 

which is obviously double counting the valence quark contribution to the first moment sum 
rule. From the practical viewpoint, this double counting is not so serious, since a^ai term 
turns out to be order of magnitude smaller than a sea- This dominance of the Dirac-sea 
contribution to the nucleon scalar charge is one of the distinguishing features of the CQSM. 
One can say that it is connected with the unique feature of this model, which is able to 
describe simultaneously a localized baryonic excitation together with the nontrivial QCD 
vacuum structure with nonzero quark condensate (or nonzero scalar quark density). In any 
case, we emphasize that the CQSM predicts fairly large scalar charge for the nucleon, i.e. 
~ 11.8. Using the current quark mass of mo 5 MeV, as an estimate, this gives 

E^iv = mo a- ~ 60 MeV, (68) 

which seems to favor relatively large values obtained from the recent analysis of the pion- 
nucleon scattering amplitude [39]-[43]. 

Next we turn to the discussion of the second moment sum rule. We first point out that 
the S{x) term in e(^=°^(x) does not contribute to the second moment. In the CQSM, then, 
the second moment of e^'^^^\x) receives contributions from two terms in the distribution, 
i.e. the valence quark term e^^^^\x) and the regular part of the vacuum polarization term 
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e^g^\x). After performing the numerical integration, we find that 

xe!ifj°\x)dx ~ 0.23, (69) 



"reg 



xe'^^=^\x)dx = y ^a;e(^=°)(a;)dx~ -0.05. (70) 

The total second moment is therefore given by 

J^^ xe^'^="\x)dx ~ 0.23 - 0.05 ~ 0.18. (71) 

We recall that, within the CQSM, there is another independent method for evaluating the 
second moment. Since we are working in the chiral limit, we rewrite (34), by setting mo = 0, 
as 

[\e^^='\x)dx^N,.^(3, (72) 
J-l Mjv 

or 



°\x)dx = Ncj^lS^ai, (73) 

f_yi='\x)dx = N^^Psea, (74) 



with 



(^vai = {Q\\{U + W)\Q), (75) 
Psea = E(n|^(C/ + C/t)|n). (76) 

These quantities id^ai ^-nd (igea can be calculated directly within the model, without invoking 
the corresponding distribution functions. Numerically, we find that 

M 

N,—-^ai - 0.23, (77) 

Nc^Psea - -0.06. (78) 

These two numbers arc consistent with the corresponding numbers in (69) and (70), obtained 
through the distribution functions, A small discrepancy between the numbers in (70) and 
(78) may be interpreted as giving a measure of numerical errors introduced by the very 
delicate interpolation method for obtaining the vacuum polarization term of e^'^^'^\x). At 
any rate, we find that the CQSM predicts relatively small but nonzero value for the second 
moment of e^'^^'^\x). Since we are working in the chiral limit (mo = 0), this appears to 
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contradict the corresponding sum rule (27) derived on the basis of the QCD equations of 
motion, which states that the second moment of e^'^^^\x) vanishes in the chiral hmit. Does 
this discrepancy simply mean the limitation of the CQSM as an effective theory of QCD? 
In our opinion, this is not necessarily the case by the following reasons. First of all, we 
point out that moment sum rules containing quark masses are somewhat delicate one, since 
the masses are generally dependent on the renormalization scale. Secondly, if the QCD 
vacuum breaks the chiral symmetry spontaneously as is generally believed, a quark acquires 
a dynamical mass of several hundred MeV, which means that massless quarks are nowhere. 
Naturally, the situation is not so simple because of the color confinement. For instance, 
according to the picture of the MIT bag model, which realizes quark confinement by hand, 
at least the vacuum inside the bag is perturbativc and the quarks inside it remains massless. 
According to Shuryak [44], the bag model is based on the idea that the hadron is a piece 
of a qualitatively different (or "perturbative" ) phase of the QCD vacuum. The physical 
picture of the CQSM for the baryon and the QCD vacuum is fairly different from that of 
the bag model. According to the words of Shuryak again, the chiral models (including the 
CQSM) assume that the vacuum is only slightly modified inside the hadron : the relative 
orientation of the right- and the left-handed quark fields is somewhat different. This last 
statement denotes the fact that, in the basic lagrangian of the CQSM, the dynamical quark 
mass parameter M appears as a product with the chiral field U^^{x), which is space-time 
dependent. It is also the cause of the fact that the product of M and P enters the r.h.s. of the 
second moment sum rule (34). This supports Schweitzer's viewpoint [20] that the quantity 
PM can be interpreted as an effective mass of quarks bound in the sohton background at 
least in the 2nd moment sum rule of e^^^^\x). Numerically, we have 

(3M ~ 51 MeV. (79) 

This value is smaller than the one obtained in [20], since the contribution of the Dirac-sea 
quarks neglected in [20] works to reduce the value of p. 

In any case, the nonzero value of the second moment of e^'^^^\x) is nothing contradictory 
at least within the framework of the CQSM in which massless quarks are nowhere at the 
model energy scale of about 600 MeV. However, we anticipate that the dynamical quark 
mass M is generally a scale dependent quantity which approaches zero in the high energy 
hmit. The naive QCD sum rule for the second moment of e^'^^^\x) would be recovered in 
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this limit. To verify the vahdity of this idea, what is crucial is experimental determination of 
the second moment sum rule at the relatively low energy scale close to the above-mentioned 
model energy scale. This may be in principle possible by inversely evoluting high energy 
data to low energy scale. 

B. isovector distribution e^'^~^\x) 

In the case of isovector distribution e^^^^-*(x), no ultraviolet regularization is needed 
because its first moment, (50), is related to the imaginary part of the Euclidian effective 
meson action in the background sohton field [45] and it is ultraviolet finite. We have checked 
that the energy level sum (43) is stable enough against the increase of the cutoff momentum 
kmax, above 12M. The final result for the isovector distribution e^'^^^\x) is shown in Fig. 6. 
The dashed curve represents the contribution of the Nc valence level quarks, the dot-dashed 
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FIG. 6: The theoretical predictions of the CQSM for e^'^^^\x). The dot-dashed curve stands for 
the contribution of Nc valence level quarks, the dashed curve the contribution of the Dirac-sea 
quarks, while the solid curves represents their sum. 

curve represents the contribution of the Dirac-sea quarks, while the solid curve represents 
their sum. In contrast to the isosinglet distribution, the Dirac-sea contribution has no 
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singularity at x — and it is a smooth function in the whole region of x. The total 
contribution is given by the solid curve. 
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FIG. 7: The comparison of the theoretical predictions for e^^^^\x) and e^^^^\x) at the model 
energy scale. 

The first moment or the x-integral of this total contribution gives the following value 

y\(^=^^(x)dx ~ 0.28, (80) 

which is order of magnitude consistent with the estimate obtained from the analysis of the 
non-electromagnetic proton-neutron mass difference. Shown in Fig. 7 are the comparison of 
our final theoretical predictions for e^^=°^ (x) and e^^=^^ (x) . One confirms that the magnitude 
of e^^^^\x) is much smaller than that of e^-^^^^x) in conformity with the large Nc relation 
(11). Combining these two distributions, we can now give final theoretical predictions for 
the chiral-odd twist-3 distribution function e"(x) of each flavor a. Shown in Fig. 8 (a) are the 
distributions for the it-quark and the u-quark, while Fig. 8(b) gives the distributions for the 
d-quark and the d-quark. 
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FIG. 8: The theoretical predictions for e"(a;), e°'(x), e"(a;), and e'^(x) at the model energy scale. 
C. comparison with the empirical information from the CLAS measurements 



Here, we make a very preliminary comparison of our theoretical predictions for e{x) with 
the empirical information extracted from the high-energy semi-inclusive scatterings. Because 
of its chiral-odd nature, the distribution function e{x) does not appear in inclusive DIS cross 
sections. To extract any information for it, we must therefore carry out more specific semi- 
inclusive type scattering experiments. Very recently, such an experiment has in fact been 
done by the CLAS Collaboration [23]. They measured the azimuthal asymmetry Alu in 
the electro-production of pions from deeply inelastic scatterings of longitudinally polarized 
electrons off unpolarized protons. 

The first theoretical analysis of the CLAS data was carried out by Efremov et al. [29], [30]. 
Their analysis assume that the beam single spin asymmetry measured by the CLAS group 
is dominantly generated by the so-called Collins mechanism [46]. Under this assumption 
together with a particular parameterization for the Collins fragmentation function, they 
were able to extract the first information on the chiral-odd twist-3 distribution function 
e{x). Recently, this analysis was criticized by Yuan [47]. According to him, there may be 
another mechanism which competes with the Collins mechanism [48], [49]. It is the leading 
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order transverse momentum dependent parton distribution hi{x, k±) convoluted with chiral- 
odd fragmentation function e{z). After all, the fact is that we still have poor knowledge 
about the mechanism that generates the beam single spin asymmetry in semi-inclusive deep 
inelastic scatterings. We must understand the mechanism of parton fragmentation processes 
into hadrons, especially the physics of time- reversal-odd fragmentation functions [46], [50]. 
We must also clarify the dynamics of transverse momentum dependent parton distribution 
functions in combination with the physics of chiral-odd fragmentation functions [48]-[50]. 
A truly reliable extraction of the chiral-odd twist-3 distribution function e{x), which is of 
our primary concern here, can be achieved only after more complete understanding of the 
above-mentioned mechanisms of semi-inclusive DIS processes. 

Keeping this fact in mind, we shall proceed here by assuming the dominance of the Collins 
mechanism. Under this assumption, the asymmetry measured by the CLAS experiment is 
interpreted to be proportional to 



with y = {P ■ q)/{P ■ I), z = {P ■ ph)/{P ■ q) and s is the invariant mass squared of the 
photon-hadron system in the notation of [29] . Ag denotes the beam helicity. The chiral- and 
T-odd twist-2 "CoUins" fragmentation function H^"-{z) gives the probability of a spinless 
or unpolarized hadron to be created from a transversely polarized scattered quark. Using 
informations on H^"'[z) from HERMES data [24], [25], one can then get direct information on 
the distribution function e{x) [29], [30]. In the CLAS experiment, the azimuthal asymmetries 
A^l^lf for the process ep — * e'n'^X were measured at ~ 1.5 GeV^. Under the dominant- 
flavor-only approximation for the fragmentation functions, the semi-inclusive tt"*" production 
measures the following combination of the distributions. 



In Fig. 9, wc make a comparison between the predictions of the CQSM for the above com- 
binations of the distributions and the corresponding empirical information extracted from 
the CLAS data by Efremov et al. [29], [30] under the assumption of Collins mechanism dom- 
inance. The theoretical distribution here corresponds to the energy scale of = 1.5 GeV^. 
The scale dependence of the distribution is taken into account by solving the leading-order 
DGLAP type equation obtained in the large Nc hmit [35]. (The starting energy scale of 




(81) 




(82) 
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this evolution is taken to be Qf^j ~ 0.30 GeV^.) The distribution e"(x) + \e'^(x) extracted 
from the CLAS data contains large errors mainly due to the large uncertainties of H^{z) 
from the HERMES data [24], [25]. Still, it was emphasized in [29], [30] that the extracted 
distribution is definitely larger than the "twist-3 bound" and about two times smaller than 
the corresponding unpolarized distribution /"(.x) at the same energy scale. One sees that 
our theoretical prediction for e^{x) + \e'^{x) is in fairly good agreement with the extracted 
behavior from the CLAS data. The relatively small magnitude of the extracted e{x) indi- 
cates that there must be a significant contribution to the nN sigma-term sum rule from 
the small x region. Whether this is due to the indicated 5-function singularity in e{x) or 
it is due to the yct-unresolved Regge behavior in the small x region is difiicult to judge at 
the present stage of study. It is highly desirable to extend the region of measurements to 
smaller x region. This is important, because the unambiguous establishment of the violation 
of the t:N sigma-term sum rule would indirectly prove the existence of a novel 5-function 
singularity in the distribution function e{x) of the nucleon, which in turn may be interpteted 
as a manifestation of the nontrivial structure of QCD vacuum in an observable of a localized 
QCD excitation, i.e. the nucleon. 

Finally, we want to make some comments on the prediction for e{x) based on the MIT 
bag model. As mentioned in [30], the bag model prediction of [28] evolved to the com- 
parable energy scale of = 1 GeV^ is in quahtative agreement with the extracted e{x) 
from the CLAS data in [30]. In our opinion, this agreement should be taken as fortuitous 
by the following reason. First, as already pointed out, the isosinglet scalar charge of the 
nucleon predicted by the MIT bag model is only about 15% of the value expected from the 
phenomenological knowledge of the tiN sigma-term. The fact is that the nucleon isoscalar 
charge is a quantity of order 1 (or order N^i more precisely) in the MIT bag model or in 
any other models which contains three valence-quark degrees of freedom only. The situation 
is totally different in the CQSM. Although the contribution of the Nc valance level quarks 
is of the same order as that of the MIT bag model, the vacuum polarization effect or the 
contribution of the Dirac-sea quarks give nearly seven times larger contribution as compared 
with that of the valence quarks, thereby reproducing the correct magnitude of the nucleon 
scalar charge or the ttN sigma-term. Unfortunately, this crucial difference between the two 
models is not refiected in the observable distribution function e{x). Since the Dirac-sea 
contribution in the CQSM is nearly saturated by the 5-function singularity, it happens that 
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FIG. 9: The theoretical prediction for e(x) = e"(x) + \e (x) in comparison with the corresponding 
empirical information extracted from the CLAS data at (Q^) = l.SGeV^ under the assumption of 
Collins mechanism dominance. 

the distributions e{x) at x 7^ predicted by the two models are not extremely different with 
each other. This is the reason why the naive MIT bag model, which fails to explain the 
magnitude of the ttN sigma-term, can reproduce the empirical distribution e{x) extracted 
from the CLAS data at least qualitatively. 

Still, we will show that there are some qualitative and observable differences between the 
predictions of the CQSM and the MIT bag model. The key observation here is that, for the 
spin-independent chiral-odd twist-3 distribution functions, the MIT bag model predicts no 
flavor dependence. That is, within the framework of the naive MIT bag model, we have 

e"(x) = e'^{x), e"(x) = e^x), (83) 

or more specifically 

e"(x) + ^e'^ix) = e'^ix) + ^e"(x). (84) 

Such equalities can be expected to hold only in the fictitious limit of 00. As is in 

fact the case with the CQSM, for a finite value of Nc, the isovector distribution e^'^^^^ (x) — 
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FIG. 10: The predictions of the CQSM for e"(x) + (l/4)e^"(x) and e^(x) + (l/4)e"(x) evolved to the 
energy scale ~ 1.5 GeV^ of the CLAS data from the initial scale of the model Ql^^ ^ 0.30 GeV^. 
Also shown for comparison is the prediction of the MIT bag model evolved to the same scale from 
somewhat lower energy scale of — 0.16 GeV^. 

e"(x) — e^{x) does not vanish, so that we definitely expect that 

e'^ix) + ^e'^ix) ^ e'^ix) + ^e^{x). (85) 

Fig. 10 shows the comparison of the predictions of the two models for the distributions 
e"(x) + |e^(a;) and e'^{x) + |e"(,T) evolved to the energy scale of CLAS experiment, i.e. 

~ 1.5 GeV^ from the initial energy scale of the model Qf^j — 0.30 GeV^. The solid and 
dashed curves here stand for the predictions of the CQSM, respectively for e^{x) + |e'^(a;) 
and e^{x) + |e"(x). On the other hand, the dot-dashed curve represents the prediction of the 
MIT bag model, which gives an identical answer for both these combinations of distributions. 
One sees that the CQSM predicts a sizably large difference between the two distributions 
e^{x) + \e'^{x) and e'^{x) + |e"(x), in sharp contrast to the MIT bag model. In principle, the 
possible differences of these two distributions can be detected by performing a comparative 
analysis of the semi-inclusive tt"*" and tt" productions. 



31 



IV. SUMMARY AND CONCLUSION 



In summary, we have given theoretical predictions for the chiral-odd twist-3 distribution 
function e"(x) of the nucleon with each flavor a on the basis of the chiral quark soliton 
model. A prominent feature of the isosinglet combination of the distributions, e^{x) + e'^{x), 
is that its first moment is proportional to the familiar ttN sigma-term and that it contains a 
delta- function singularity at a; = 0. In the previous study based on the derivative expansion 
technique, we demonstrated that the physical origin of this singularity can be traced back to 
the long-range quark-quark correlation of scalar type, which signals the spontaneous chiral 
symmetry breaking of the QCD vacuum. The present calculation, without recourse to the 
derivative expansion type approximation, has revealed the following facts. The isosinglet 
distribution e"(a;) + e^{x) consists of two parts, i.e. the contribution of Nc valence level 
quarks and that of the Dirac sea quarks in the hedgehog mean field. The former takes a 
familiar shape of distribution which has a peak around the value of x ~ 1/3. On the other 
hand, the latter certainly contains a (^-function type singularity at x — 0, but it also has 
nontrivial support for x 7^ 0. The isovector distribution e"(x) — e^{x) also consists of the 
valence and Dirac-sea contributions. For this distribution, however, no delta-function type 
singularity is observed, which means that it is a regular function in all the range of x. 

The moment sum rules of e{x) provide us with valuable information concerning the basic 
dynamical content of the model in view of the underlying theory, i.e. QCD. We showed 
that the first moment sum rule for e"{x) + e'^{x) is satisfied within the model, if and only if 
the delta-function singularity is properly taken into account. Note however that the delta- 
function term alone does not saturate the first moment or the ttN sigma-term sum rule in 
contrast to the previous argument based on the framework of the perturbative QCD. We 
also pointed out that the second moment sum rule for e^{x) + e'^{x) does not vanish even in 
the chiral limit in contrast to the QCD-equation-of-motion argument. In our opinion, this 
violation of the second moment sum rule does not necessarily show a defect of the model. It 
is rather to be interpreted as showing the limitation of the perturbative analysis as a tool of 
handling a bound state problem and/or the problem of masses nonperturbatively generated 
by the mechanism of the spontaneous chiral-symmetry breaking. We have also shown that 
the model prediction for the first moment of the isovector distribution e"(a;) — e'^{x) comes 
out to be order of magnitude consistent with the phenomenological estimate obtained from 
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the nonelectromagnetic neutron-proton mass difference. 

It was sfiown tliat tfie tfieoretical predictions for the distribution e"(x) + |e''(x) are in a 
good agreement with the corresponding empirical information extracted from the CLAS data 
for the semi-inclusive 7r+ production under the assumption of the Collins mechanism domi- 
nance. This agreement, combined with our analysis explained in the text, implies the exis- 
tence of 5-function singularity at x = in the isosinglet distribution e^{x) + e'^{x), although 
a definite conclusion must awaits for more complete measurements and more thorough un- 
derstanding of the reaction mechanism that generates the beam single spin asymmetry in 
the semi-inclusive pion productions. 

Finally, we compare our theoretical predictions with those of the MIT bag model. As 
shown in the body of the paper, the two models give accidentally close predictions for 
the distribution function e"(x) -|- \e'^{x) at x ^ 0. We have shown, however, that the 
CQSM predicts a sizably large difference between the two distributions e"(x) -|- |e''(x) and 
e'''(x) + |e"(x), for which the MIT bag model makes no difference. The predicted sizable 
difference between the two combinations of distributions will be detected by performing a 
comparative experimental analysis of the semi-inclusive tt"*" and 7r~ productions. 
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